System and method for automated seismic interpretation

ABSTRACT

A computer-implemented method is described for automated seismic interpretation that includes receiving, at one or more processors, a pre-stack seismic dataset representative of the subsurface volume of interest; performing, via the one or more processors, a machine learning algorithm on the pre-stack seismic dataset to identify seismic facies based on AVA clusters, wherein the identified seismic facies include cluster sequences in depth for a plurality of spatial x-y locations in the subsurface volume; performing, via the one or more processors, seismic interpretation of the seismic dataset based on the identified seismic facies to generate a digital image of the seismic interpretation; and displaying the digital image of the seismic interpretation on a user interface.

CROSS-REFERENCE TO RELATED APPLICATIONS

The present application is a continuation-in-part of U.S. application Ser. No. 15/730,814, filed Oct. 12, 2017, which claims the priority benefit of U.S. Provisional Application Ser. No. 62/408,188, filed on Oct. 14, 2016, both of which are incorporated herein by reference in their entirety.

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT

Not applicable.

TECHNICAL FIELD

The disclosed embodiments relate generally to techniques for automated seismic interpretation of subsurface volumes of interest and, in particular, to a method of identifying seismic facies and using them for automated interpretation of a subsurface volume of interest including complex geology.

BACKGROUND

Seismic exploration involves surveying subterranean geological media for hydrocarbon deposits. A survey typically involves deploying seismic sources and seismic sensors at predetermined locations. The sources generate seismic waves, which propagate into the geological medium creating pressure changes and vibrations. Variations in physical properties of the geological medium give rise to changes in certain properties of the seismic waves, such as their direction of propagation and other properties.

Portions of the seismic waves reach the seismic sensors. Some seismic sensors are sensitive to pressure changes (e.g., hydrophones), others to particle motion (e.g., geophones), and industrial surveys may deploy one type of sensor or both. In response to the detected seismic waves, the sensors generate corresponding electrical signals, known as traces, and record them in storage media as seismic data. Seismic data will include a plurality of “shots” (individual instances of the seismic source being activated), each of which are associated with a plurality of traces recorded at the plurality of sensors.

Seismic data is processed to create seismic images that can be interpreted to identify subsurface geologic features including hydrocarbon deposits. The seismic image can be used to obtain seismic facies. Seismic facies are groups of reflections in the seismic image that can be categorized based on characteristics such as amplitude (e.g., amplitude variation with offset/angle), continuity, geometry, and/or texture. Knowledgeable practitioners can interpret these facies as corresponding to subsurface properties such as lithology, depositional environment, and fluid content. This information may be used in a hydrocarbon exploration setting to search for ideal geological settings for hydrocarbon reservoirs, and for model-building purposes to assign appropriate properties (velocity, density, permeability, porosity, etc.) within the individual facies.

Prior art includes approaches based on classifying data on a trace-by-trace basis using, for example, neural nets or self-organizing maps. This can be detrimental to accurately representing the true geometry or connectivity of geological features. Moreover, these trace-by-trace methods are time consuming and create bottlenecks in the process of characterizing the subsurface.

The ability to identify facies in the subsurface is crucial to our ability to make the most appropriate choices for purchasing materials, operating safely, and successfully completing projects. Project cost is dependent upon accurate prediction of the position of physical boundaries within the Earth. Decisions include, but are not limited to, budgetary planning, obtaining mineral and lease rights, signing well commitments, permitting rig locations, designing well paths and drilling strategy, preventing subsurface integrity issues by planning proper casing and cementation strategies, and selecting and purchasing appropriate completion and production equipment.

There exists a need for determining 3-D distributions of seismic facies from seismic images more efficiently and accurately that will allow better seismic interpretation of potential hydrocarbon reservoirs.

SUMMARY

In accordance with some embodiments, a method of automated seismic interpretation that includes receiving, at one or more processors, a pre-stack seismic dataset representative of the subsurface volume of interest; performing, via the one or more processors, a machine learning algorithm on the pre-stack seismic dataset to identify seismic facies based on AVA clusters, wherein the identified seismic facies include cluster sequences in depth for a plurality of spatial x-y locations in the subsurface volume; performing, via the one or more processors, seismic interpretation of the seismic dataset based on the identified seismic facies to generate a digital image of the seismic interpretation; and displaying the digital image of the seismic interpretation on a user interface is disclosed.

In another aspect of the present invention, to address the aforementioned problems, some embodiments provide a non-transitory computer readable storage medium storing one or more programs. The one or more programs comprise instructions, which when executed by a computer system with one or more processors and memory, cause the computer system to perform any of the methods provided herein.

In yet another aspect of the present invention, to address the aforementioned problems, some embodiments provide a computer system. The computer system includes one or more processors, memory, and one or more programs. The one or more programs are stored in memory and configured to be executed by the one or more processors. The one or more programs include an operating system and instructions that when executed by the one or more processors cause the computer system to perform any of the methods provided herein.

BRIEF DESCRIPTION OF THE DRAWINGS

The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the Office upon request and payment of the necessary fee.

FIG. 1 illustrates a flowchart of a method of seismic facies identification, in accordance with some embodiments;

FIG. 2 illustrates a flowchart of a method of seismic facies identification, in accordance with other embodiments;

FIG. 3 illustrates a flowchart of a method of seismic facies identification, in accordance with other embodiments;

FIG. 4 illustrates a flowchart of a method of seismic facies identification, in accordance with other embodiments;

FIG. 5 demonstrates the result of one of the embodiments;

FIGS. 6A-6D demonstrate detailed steps of one of the embodiments;

FIG. 7 is a block diagram illustrating a seismic facies identification system, in accordance with some embodiments;

FIG. 8 illustrates a flowchart of a method of automated seismic interpretation, in accordance with some embodiments;

FIGS. 9A-9C demonstrate steps and results of an embodiment;

FIG. 10 illustrates steps of an embodiment; and

FIG. 11 illustrates steps of an embodiment.

Like reference numerals refer to corresponding parts throughout the drawings.

DETAILED DESCRIPTION OF EMBODIMENTS

Described below are methods, systems, and computer readable storage media that provide a manner of seismic facies identification via machine learning. These embodiments are designed to be of particular use for seismic imaging of subsurface volumes in geologically complex areas. These approaches use geometric and/or textural characteristics associated with seismic facies and identifies them using established or novel machine learning techniques. Machine learning techniques may be able to take advantage of significant computational processing power to uncover patterns in data that would be difficult for human interpreters to spot, leading to more accurate facies classifications.

Reference will now be made in detail to various embodiments, examples of which are illustrated in the accompanying drawings. In the following detailed description, numerous specific details are set forth in order to provide a thorough understanding of the present disclosure and the embodiments described herein. However, embodiments described herein may be practiced without these specific details. In other instances, well-known methods, procedures, components, and mechanical apparatus have not been described in detail so as not to unnecessarily obscure aspects of the embodiments.

The present invention includes embodiments of a method and system for automated seismic interpretation based on seismic facies identification that produces 3-D digital seismic images. This will help to alleviate bottlenecks created by the necessity of manually interpreting seismic facies for exploration, model-building, or reservoir modeling purposes. The present invention may be able to take advantage of significant computational processing power uncover patterns in data that would be difficult for human interpreters to spot, leading to more accurate seismic interpretation.

Advantageously, those of ordinary skill in the art will appreciate, for example, that the embodiments provided herein may be utilized to generate a more accurate digital seismic image (i.e., the corrected digital seismic image). The more accurate digital seismic image may improve hydrocarbon exploration and improve hydrocarbon production. The more accurate digital seismic image may provide details of the subsurface that were illustrated poorly or not at all in traditional seismic images. Moreover, the more accurate digital seismic image may better delineate where different features begin, end, or any combination thereof. As one example, the more accurate digital seismic image may illustrate faults more accurately. As another example, assume that the more accurate digital seismic image indicates the presence of a hydrocarbon deposit. The more accurate digital seismic image may delineate more accurately the bounds of the hydrocarbon deposit so that the hydrocarbon deposit may be produced.

Those of ordinary skill in the art will appreciate, for example, that the more accurate digital seismic image may be utilized in hydrocarbon exploration and hydrocarbon production for decision making. For example, the more accurate digital seismic image may be utilized to pick a location for a wellbore. Those of ordinary skill in the art will appreciate that decisions about (a) where to drill one or more wellbores to produce the hydrocarbon deposit, (b) how many wellbores to drill to produce the hydrocarbon deposit, etc. may be made based on the more accurate digital seismic image. The more accurate digital seismic image may even be utilized to select the trajectory of each wellbore to be drilled. Moreover, if the delineation indicates a large hydrocarbon deposit, then a higher number of wellbore locations may be selected and that higher number of wellbores may be drilled, as compared to delineation indicating a smaller hydrocarbon deposit.

Those of ordinary skill in the art will appreciate, for example, that the more accurate digital seismic image may be utilized in hydrocarbon exploration and hydrocarbon production for control. For example, the more accurate digital seismic image may be utilized to steer a tool (e.g., drilling tool) to drill a wellbore. A drilling tool may be steered to drill one or more wellbores to produce the hydrocarbon deposit. Steering the tool may include drilling around or avoiding certain subsurface features (e.g., faults, salt diapirs, shale diapirs, shale ridges, pockmarks, buried channels, gas chimneys, shallow gas pockets, and slumps), drilling through certain subsurface features (e.g., hydrocarbon deposit), or any combination thereof depending on the desired outcome. As another example, the more accurate digital seismic image may be utilized for controlling flow of fluids injected into or received from the subsurface, the wellbore, or any combination thereof. As another example, the more accurate digital seismic image may be utilized for controlling flow of fluids injected into or received from at least one hydrocarbon producing zone of the subsurface. Chokes or well control devices, positioned on the surface or downhole, may be used to control the flow of fluid into and out. For example, certain subsurface features in the more accurate digital seismic image may prompt activation, deactivation, modification, or any combination thereof of the chokes or well control devices so as control the flow of fluid. Thus, the more accurate digital seismic image may be utilized to control injection rates, production rates, or any combination thereof.

Those of ordinary skill in the art will appreciate, for example, that the more accurate digital seismic image may be utilized to select completions, components, fluids, etc. for a wellbore. A variety of casing, tubing, packers, heaters, sand screens, gravel packs, items for fines migration, etc. may be selected for each wellbore to be drilled based on the more accurate digital seismic image. Furthermore, one or more recovery techniques to produce the hydrocarbon deposit may be selected based on the more accurate digital seismic image.

In short, those of ordinary skill in the art will appreciate that there are many decisions (e.g., in the context of (a) steering decisions, (b) landing decisions, (c) completion decisions, (d) engineering control systems and reservoir monitoring in the following but not limited to: Tow Streamer, Ocean Bottom Sensor, VSP, DASVSP, and imaging with both primaries and free surface multiple, etc.) to make in the hydrocarbon industry and making proper decisions based on more accurate digital seismic images should improve the likelihood of safe and reliable operations. For simplicity, the many possibilities, including wellbore location, component selection for the wellbore, recovery technique selection, controlling flow of fluid, etc., may be collectively referred to as managing a subsurface reservoir.

FIG. 1 illustrates a flowchart of a method 100 for seismic facies identification in a subsurface volume of interest. At operation 10, a seismic dataset is received. As previously described, the seismic dataset includes a plurality of traces recorded at a plurality of seismic sensors. This dataset may have already been subjected to a number of seismic processing steps, such as deghosting, multiple removal, spectral shaping, imaging, and the like. These examples are not meant to be limiting. Those of skill in the art will appreciate that there are a number of useful seismic processing steps that may be applied to seismic data. The seismic dataset can be raw trace data, pre-stack gathers, or a migrated image/volume that may be pre-stack or post-stack; in the case of migrated data, the input may or may not be volumetrically flattened prior to analysis. In the case of pre-stack gathers or a pre-stack migrated image, the dataset will include an offset and/or angle dimension, representative of the distance between the source-receiver pairs. The offset and/or angle dimension will include seismic events with amplitudes that vary along the offset or angle which those of skill in the art know as amplitude variation with offset (AVO) or amplitude variation with angle (AVA). For the purposes of this paper, AVO and AVA are collectively referred to as AVA.

Continuing method 100, at operation 12 the seismic dataset is subjected to unsupervised machine learning. The input data is automatically classified or clustered into distinct groups based on characteristics of the input data. This classification can be accomplished by algorithms such as, but not limited to, k-means clustering (described in Hartigan and Wong, 1979), self-organizing maps (Kohonen, 1990), or deep learning techniques (Lecun et al., 2015). The classification allows the input data to be aggregated into groups that improve the signal-to-noise ratio in appropriate ways for those clusters. In an embodiment, the unsupervised classification may be a data-driven, automated approach to seismic aggregation (i.e., clustering) capable of automatically determining the appropriate data to evaluate for robust, representative AVA responses corresponding to a subsurface location. The classification may be performed on multiple local and global scales to determine regions of rock and fluid property similarity and differences. For example, local clustering may aggregate localized data based on raw amplitudes or seismic-derived features such as texture attributes. Global clustering may be done for refinement of clustering for specific AVA Class response or additional rock and fluid property interrogation (e.g., spatial porosity variations, volumetric prospect screening), possibly distinguishing between different fluid content (e.g., brine vs. gas) in the same lithology. In another embodiment, unsupervised machine learning can cluster seismic data with a consistent, stable AVA Class response based on the intercept, gradient, and polarity of the seismic events. An example of clustering based on the AVA responses is demonstrated in FIGS. 6A-6D.

FIG. 6A demonstrates the input data used for the clustering. In FIG. 6A, panel 60 shows map view depth slices through a seismic image volume for a near stack, mid stack, and far stack image. Those of skill in the art will understand that the partial stacks (near, mid, and far) are a summation of traces along the offset or angle axis for a range of offsets/angles representative of source-receiver distances that are small (near), mid-range (mid), and large (far), calculated from a pre-stack migration. The ten panels labeled “Local AVA Response” in FIG. 6A show examples of amplitude variation with angle (AVA) for individual sample points in the images. Locally, at any subsurface sample the AVA response can be noisy and there is much uncertainty associated with this one measurement.

FIG. 6B demonstrates “soft” clustering of the various local AVA responses, aggregating the local AVA responses into soft clusters 61, 62, 63, and 64 based on the type of AVA response. In this soft clustering, the clusters map overlap somewhat and the individual sample points may belong to more than one cluster. Soft clustering may be accomplished, by way of example and not limitation, by Fuzzy C-Means Clustering.

FIG. 6C demonstrates “hard” clustering of the various local AVA responses, aggregating the local AVA responses into clusters 65, 66, 67, and 68 based on the type of AVA response and the probability of it belonging to a cluster. In hard clustering, each sample will belong to just one cluster. Once the local AVA responses are clustered, a global statistical AVA response for each cluster can be calculated from the local AVA responses in that cluster.

FIG. 6D panel 70 shows the map view images of the samples divided into clusters for three horizons. FIG. 6D panel 71 shows the map view images of the probability of the samples for belonging to one particular cluster for three horizons. One of skill in the art will recognize that these map view images display a recognizable geologic feature known as a meandering stream. Such recognizable geologic features confirm that the clustering is likely to have been accurate.

After the unsupervised machine learning has generated the distinct groups, the groups may optionally be labeled as various facies at operation 14. Whether or not the facies are labeled, the distinct groups are then identified as 3-D geobodies, meaning that the geometry of each group is defined in extent and volume, at operation 16. If desired, the method may repeat (i.e. iterate) operations 12, 14, and 16 until a user-specified stopping point, such as a specified number of iterations or a 3-D volume wherein all points are included within a geobody.

An alternate embodiment is illustrated as method 120 in FIG. 2. In this embodiment, a seismic dataset is received at operation 20, much like operation 10 of method 100. An optional additional input of a set of facies labels may also be received at operation 22. The facies labels consist of a sub-volume of the input data that has been wholly or partly interpreted and manually assigned a seismic facies category. At the locations where such interpretation has been performed, the input data will therefore consist of both the seismic amplitude at each pixel or voxel of the image, and the corresponding facies label at that pixel or voxel.

The seismic dataset and optional set of facies labels are input to supervised machine learning at operation 24. In this case, the labeled training dataset is first used to “train” a classifier; the classifier is then applied to the unlabeled input data in order to produce a facies distribution that is now specifically labeled. Classification algorithms in this case can include decision-tree methods such as Random Forest (Breiman, 2001); classical neural net methods (Lippmann, 1987); or more novel deep learning methods such as a convolutional neural net (Krizhevsky et al., 2012). Method 120 may optionally include using artificial intelligence computations in the data space to recognize stratigraphic patterns such as channel features. Training data for this optional step may be obtained from sources such as real or synthetic well logs, or simulated stratigraphic models. The 3-D geometry of the geobodies is defined at operation 28. This embodiment may also iterate over operations 24, 26, and 28.

A specific example of method 120 is shown in FIG. 4 as method 160. It receives a seismic image at operation 40 and flattens it at operation 42, using, for example, a method described in U.S. Pat. No. 7,769,546, Method for Indexing a Subsurface Volume For The Purpose of Inferring Geologic Information, or U.S. patent application Ser. No. 14/595,964, System and Method for Generating a Depositional Sequence Volume from Seismic Data. The flattened image is then interpreted and used as a partial image at operation 44 for input to the Random Forest algorithm at operation 48 to build a classifier based on limited labels provided by an interpreter. The classifier can then be applied to a full seismic volume to produce a 3-D image with the geobodies identified.

Method 160 optionally includes calculating textural attributes at operation 46 as they in particular may be useful for defining the facies distributions geometrically, as described above. The textural attribute calculations result in multiple attribute images that could be input to the Random Forest algorithm at operation 48 to produce a 3-D image with the geobodies identified.

FIG. 3 shows an optional workflow extension (method 140) of any of the methods 100, 120, or 160. If an appropriate computational stratigraphy model is available, it can be received as input 30 to enhance the facies distribution result by inferring additional information about the lithology, connectivity, and/or flow characteristics of certain facies at operation 32. This can be achieved via an additional matching step at operation 34, in which a machine learning technique matches the seismic expression of specific features from the model to patterns found in the input seismic data. This extension is viable for either the supervised or unsupervised learning cases discussed above.

FIG. 5 demonstrates the result of using one of the embodiments of the present invention. In particular, it shows the result of method 120 of FIG. 2. Panel 50 shows a single depth slice of a seismic image. The input facies labels are shown as the red line 51, green rectangle 53, and purple rectangle and line 54. After performing the method 120 of FIG. 2, this depth slice now shows the classification result 52 wherein the seismic facies have been identified as extensive geobodies.

The methods described in FIGS. 1 through 4 may be extended to enable automated seismic interpretation as shown in FIG. 8. In conventional seismic interpretation, a human interpreter picks seismic horizons, including attempting to match horizons that may be discontinuous due to faults, pinch-outs, unconformities, and the like. Since humans cannot identify seismic facies as presented in the embodiments shown in methods 100, 120, 140, and 160, the interpretation may be inaccurate especially in areas of complex geology. The present embodiment of method 1100 overcomes these difficulties by performing a computer-implemented interpretation. In FIG. 8, method 1100 begins at operation 80 by receiving seismic data. At operation 82, the seismic facies are identified, for example using any of the methods 100, 120, 140, or 160 as previously described. However, in the present embodiment, the seismic facies are identified as clusters in vertical sequences. Each vertical sequence is constructed for a different x-y spatial location. The vertical sequences may be generated for every x-y location or at a user-specified density in x and/or y. Additionally, the vertical density (in time or depth) of the sequences may be user-specified. Based on these vertical sequences of seismic clusters, the method may move on to operation 86 to perform automated seismic interpretation based on the seismic facies. During the automated seismic interpretation, the vertical sequences of seismic clusters are used to trace the seismic horizons in the x-y domain as well as the depth (or time) domain. This is of particular use in complex geologic volumes including faults, where the seismic horizon may be vertically offset by the fault; by relying on the seismic facies clusters, a horizon with the same seismic facies can be identified across the fault, preventing erroneous interpretations that misjudge the displacement across the fault. FIGS. 9A through 10 demonstrate the method.

Method 1100 includes an optional operation 84 for fluid substitution of the seismic facies. This operation is also of particular use in geologically complex volumes, where the fluid content of a rock layer may be different on opposites sides of a fault, caused, for example and not limited to, by fluid migration along the fault or trapping at the fault. This operation is based on the assumption that there is a baseline lithology/brine response of a particular rock layer (e.g., the reservoir layer). The P50 AVA response of this baseline is associated with the clusters that have subparallel P50 AVA responses (see, for example, FIGS. 10 and 11). These subparallel clusters are either gas or oil and have the same lithology. This is the case of clusters 4 and 5 in the FIGS. 10 and 11. By substituting different fluids, it is possible for operation 86 of method 1100 to correctly match seismic facies representative of the same lithology across the faults despite the different fluids.

Once the automated seismic interpretation has completed, the results may be displayed on a user interface at operation 88. The seismic interpretation may be displayed by itself or overlain on a seismic data volume. The display may be 2-D or 3-D.

Steps of method 1100 are illustrated by FIGS. 9A-11. In FIG. 9A, a slice of a mid-stack seismic volume is shown. This is a partial stack of the mid-range angles of a seismic image. In FIG. 9B, this mid-stack seismic volume has been overlain with seismic horizons that have been automatically interpreted as described by method 1100. The horizons are labeled based on the AVA cluster to which they were assigned. The representative AVA responses of the clusters is shown in FIG. 9C.

FIG. 10 shows the same mid-stack seismic volume along with the vertical sequences of clusters, as seen on either side of Fault A. Identifying the clusters on either side of the fault allows the automated interpretation of horizons on either side of the fault.

FIG. 11 illustrates the use of operation 84 of method 1100. Here the interpretation correlates sequences across faults. Note that Cluster 4 exhibits a fluid effect indicative of hydrocarbon (subparallel AVA trend to Cluster 5 and higher amplitude up dip). Cluster logic (correlation of lithologies across faults that are the same but contain different fluids) will associate Cluster 4 with Cluster 5. The 3D reconstruction process will use these automatic correlations across faults. In this embodiment, the method not only produces an interpretation of the structure of the subsurface but also provides information on the possible fluid content, enabling the identification of potential hydrocarbon reservoirs.

FIG. 7 is a block diagram illustrating a seismic facies identification system 500, in accordance with some embodiments. While certain specific features are illustrated, those skilled in the art will appreciate from the present disclosure that various other features have not been illustrated for the sake of brevity and so as not to obscure more pertinent aspects of the embodiments disclosed herein.

To that end, the seismic facies identification system 500 includes one or more processing units (CPUs) 502, one or more network interfaces 508 and/or other communications interfaces 503, memory 506, and one or more communication buses 504 for interconnecting these and various other components. The seismic facies identification system 500 also includes a user interface 505 (e.g., a display 505-1 and an input device 505-2). The communication buses 504 may include circuitry (sometimes called a chipset) that interconnects and controls communications between system components. Memory 506 includes high-speed random access memory, such as DRAM, SRAM, DDR RAM or other random access solid state memory devices; and may include non-volatile memory, such as one or more magnetic disk storage devices, optical disk storage devices, flash memory devices, or other non-volatile solid state storage devices. Memory 506 may optionally include one or more storage devices remotely located from the CPUs 502. Memory 506, including the non-volatile and volatile memory devices within memory 506, comprises a non-transitory computer readable storage medium and may store seismic data, velocity models, seismic images, and/or geologic structure information.

In some embodiments, memory 506 or the non-transitory computer readable storage medium of memory 506 stores the following programs, modules and data structures, or a subset thereof including an operating system 516, a network communication module 518, and a seismic imaging module 520.

The operating system 516 includes procedures for handling various basic system services and for performing hardware dependent tasks.

The network communication module 518 facilitates communication with other devices via the communication network interfaces 508 (wired or wireless) and one or more communication networks, such as the Internet, other wide area networks, local area networks, metropolitan area networks, and so on.

In some embodiments, the seismic facies module 520 executes the operations of method 100. Seismic facies module 520 may include data sub-module 525, which handles the seismic dataset including seismic gathers 525-1 through 525-N. This seismic data is supplied by data sub-module 525 to other sub-modules.

Machine learning sub-module 522 contains a set of instructions 522-1 and accepts metadata and parameters 522-2 that will enable it to execute, for example, operation 12 of method 100, operation 24 of method 120, operation 32 of method 140, and/or operation 44 of method 160. The labeling function sub-module 523 contains a set of instructions 523-1 and accepts metadata and parameters 523-2 that will enable it to contribute to, for example, operation 14 of method 100, operation 22 of method 120, and/or operation 32 of method 140. The geometry sub-module 524 contains a set of instructions 524-1 and accepts metadata and parameters 524-2 that will enable it to execute, for example, operation 16 of method 100, operation 28 of method 120, operation 34 of method 140, and/or operations 46 and 48 of method 160. Although specific operations have been identified for the sub-modules discussed herein, this is not meant to be limiting. Each sub-module may be configured to execute operations identified as being a part of other sub-modules, and may contain other instructions, metadata, and parameters that allow it to execute other operations of use in processing seismic data and generate the seismic facies image. For example, any of the sub-modules may optionally be able to generate a display that would be sent to and shown on the user interface display 505-1. In addition, any of the seismic data or processed seismic data products may be transmitted via the communication interface(s) 503 or the network interface 508 and may be stored in memory 506.

Method 100 is, optionally, governed by instructions that are stored in computer memory or a non-transitory computer readable storage medium (e.g., memory 506 in FIG. 5) and are executed by one or more processors (e.g., processors 502) of one or more computer systems. The computer readable storage medium may include a magnetic or optical disk storage device, solid state storage devices such as flash memory, or other non-volatile memory device or devices. The computer readable instructions stored on the computer readable storage medium may include one or more of: source code, assembly language code, object code, or another instruction format that is interpreted by one or more processors. In various embodiments, some operations in each method may be combined and/or the order of some operations may be changed from the order shown in the figures. For ease of explanation, methods 100, 120, 140, and 160 are described as being performed by a computer system, although in some embodiments, various operations of method 100 are distributed across separate computer systems.

The system 500 may also be extended to include modules capable of executing the operations of method 1100. The modules may include a fluid substitution module and an interpretation module.

While particular embodiments are described above, it will be understood it is not intended to limit the invention to these particular embodiments. On the contrary, the invention includes alternatives, modifications and equivalents that are within the spirit and scope of the appended claims. Numerous specific details are set forth in order to provide a thorough understanding of the subject matter presented herein. But it will be apparent to one of ordinary skill in the art that the subject matter may be practiced without these specific details. In other instances, well-known methods, procedures, components, and circuits have not been described in detail so as not to unnecessarily obscure aspects of the embodiments.

The terminology used in the description of the invention herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used in the description of the invention and the appended claims, the singular forms “a,” “an,” and “the” are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will also be understood that the term “and/or” as used herein refers to and encompasses any and all possible combinations of one or more of the associated listed items. It will be further understood that the terms “includes,” “including,” “comprises,” and/or “comprising,” when used in this specification, specify the presence of stated features, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, operations, elements, components, and/or groups thereof.

As used herein, the term “if” may be construed to mean “when” or “upon” or “in response to determining” or “in accordance with a determination” or “in response to detecting,” that a stated condition precedent is true, depending on the context. Similarly, the phrase “if it is determined [that a stated condition precedent is true]” or “if [a stated condition precedent is true]” or “when [a stated condition precedent is true]” may be construed to mean “upon determining” or “in response to determining” or “in accordance with a determination” or “upon detecting” or “in response to detecting” that the stated condition precedent is true, depending on the context.

Although some of the various drawings illustrate a number of logical stages in a particular order, stages that are not order dependent may be reordered and other stages may be combined or broken out. While some reordering or other groupings are specifically mentioned, others will be obvious to those of ordinary skill in the art and so do not present an exhaustive list of alternatives. Moreover, it should be recognized that the stages could be implemented in hardware, firmware, software or any combination thereof.

The foregoing description, for purpose of explanation, has been described with reference to specific embodiments. However, the illustrative discussions above are not intended to be exhaustive or to limit the invention to the precise forms disclosed. Many modifications and variations are possible in view of the above teachings. The embodiments were chosen and described in order to best explain the principles of the invention and its practical applications, to thereby enable others skilled in the art to best utilize the invention and various embodiments with various modifications as are suited to the particular use contemplated.

REFERENCES

-   U.S. Pat. No. 7,769,546, Method for Indexing a Subsurface Volume For     The Purpose of Inferring Geologic Information -   U.S. patent application Ser. No. 14/595,964, System and Method for     Generating a Depositional Sequence Volume from Seismic Data -   Breiman, Leo. “Random forests.” Machine learning 45, no. 1 (2001):     5-32. -   Hartigan, John A., and Manchek A. Wong. “Algorithm AS 136: A k-means     clustering algorithm.” Journal of the Royal Statistical Society.     Series C (Applied Statistics) 28, no. 1 (1979): 100-108. -   Kohonen, Teuvo. “The self-organizing map.” Proceedings of the IEEE     78, no. 9 (1990): 1464-1480. -   Krizhevsky, Alex, Ilya Sutskever, and Geoffrey E. Hinton. “Imagenet     classification with deep convolutional neural networks.” In Advances     in neural information processing systems, pp. 1097-1105. 2012. -   LeCun, Yann, Yoshua Bengio, and Geoffrey Hinton. “Deep learning.”     Nature 521, no. 7553 (2015): 436-444. -   Lippmann, Richard. “An introduction to computing with neural nets.”     IEEE Assp magazine 4, no. 2 (1987): 4-22. 

What is claimed is:
 1. A computer-implemented method of automated seismic interpretation of a subsurface volume including complex geology, comprising: receiving, via a computer processor, a pre-stack seismic dataset representative of a subsurface volume of interest; identifying, via the computer processor, seismic facies in the pre-stack seismic dataset based on machine learning, wherein the machine learning identifies the seismic facies based on clustering of amplitude variation with angle (AVA) in the pre-stacked seismic dataset, wherein the seismic facies correspond to AVA clusters and include cluster sequences in depth for a plurality of spatial x-y locations in the subsurface volume of interest; generating, via the computer processor, an ensemble of fluid-related AVA responses for the AVA clusters based on a fluid substitution operation; identifying, via the computer processor, based on the ensemble of fluid-related AVA responses, related seismic facies across faults indicated in the pre-stack seismic dataset, wherein the related seismic facies contain different fluids but similar lithologies than the seismic facies; performing, via the computer processor, seismic interpretation of the pre-stacked seismic dataset based on the seismic facies and the related seismic facies to generate a digital image corresponding to the seismic interpretation; and displaying, via the computer processor, the digital image on a user interface.
 2. The computer-implemented method of claim 1, wherein the machine learning is unsupervised.
 3. The computer-implemented method of claim 2, wherein the machine learning is a unsupervised classifier that clusters the pre-stacked seismic dataset based at least on one of seismic amplitudes or seismic-derived features, wherein the seismic-derived features include at least texture attributes.
 4. The computer-implemented method of claim 1, wherein the machine learning is supervised.
 5. The computer-implemented method of claim 1, further comprising using the seismic interpretation for seismic volume reconstruction.
 6. The computer-implemented method of claim 5, wherein the seismic volume reconstruction reconstructs the subsurface volume of interest around the faults.
 7. The computer-implemented method of claim 1, wherein the machine learning identifies the seismic facies based on clustering of amplitude variation with offset (AVO) in the pre-stacked seismic dataset, and wherein the seismic facies correspond to AVO clusters.
 8. The computer-implemented method of claim 7, further comprising: generating an ensemble of fluid-related AVO responses for the AVO clusters based on a fluid substitution operation; and identifying, based on the ensemble of fluid-related AVO responses, related seismic facies across faults indicated in the pre-stack seismic dataset, wherein the related seismic facies contain different fluids but similar lithologies than the seismic facies.
 9. The computer-implemented method of claim 1, wherein the pre-stacked seismic dataset includes data relating to a plurality of traces recorded via a plurality of seismic sensors.
 10. The computer-implemented method of claim 1, wherein the data relating to the plurality of traces has been pre-processed using at least one of deghosting, multiple removal, or spectral shaping.
 11. A system, comprising: one or more processors; memory; and one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, the one or more programs including instructions that when executed by the one or more processors cause the system to: receive a pre-stack seismic dataset representative of a subsurface volume of interest; identify seismic facies in the pre-stack seismic dataset based on machine learning, wherein the machine learning identifies the seismic facies based on clustering of amplitude variation with angle (AVA) in the pre-stacked seismic dataset, wherein the seismic facies correspond to AVA clusters and include cluster sequences in depth for a plurality of spatial x-y locations in the subsurface volume of interest; generate an ensemble of fluid-related AVA responses for the AVA clusters based on a fluid substitution operation; identify, based on the ensemble of fluid-related AVA responses, related seismic facies across faults indicated in the pre-stack seismic dataset, wherein the related seismic facies contain different fluids but similar lithologies than the seismic facies; perform seismic interpretation of the pre-stacked seismic dataset based on the seismic facies and the related seismic facies to generate a digital image of the seismic interpretation; and display the digital image on a user interface.
 12. The system of claim 11, wherein the machine learning identifies the seismic facies based on clustering of amplitude variation with offset (AVO) in the pre-stacked seismic dataset, and wherein the seismic facies correspond to AVO clusters.
 13. The system of claim 12, wherein the one or more programs, when executed, further cause the system to: generate an ensemble of fluid-related AVO responses for the AVO clusters based on a fluid substitution operation; and identify, based on the ensemble of fluid-related AVO responses, related seismic facies across faults indicated in the pre-stack seismic dataset, wherein the related seismic facies contain different fluids but similar lithologies than the seismic facies.
 14. The system of claim 11, wherein the pre-stacked seismic dataset includes data relating to a plurality of traces recorded via a plurality of seismic sensors.
 15. The system of claim 11, wherein the data relating to the plurality of traces has been pre-processed using at least one of deghosting, multiple removal, or spectral shaping.
 16. A non-transitory computer readable storage medium storing one or more programs, the one or more programs comprising instructions that, when executed by one or more processors of a computing system, cause the computing system to: receive a pre-stack seismic dataset representative of a subsurface volume of interest; identify seismic facies in the pre-stack seismic dataset based on machine learning, wherein the machine learning identifies the seismic facies based on clustering of amplitude variation with angle (AVA) in the pre-stacked seismic dataset, wherein the seismic facies correspond to AVA clusters and include cluster sequences in depth for a plurality of spatial x-y locations in the subsurface volume of interest; generate an ensemble of fluid-related AVA responses for the AVA clusters based on a fluid substitution operation; identify, based on the ensemble of fluid-related AVA responses, related seismic facies across faults indicated in the pre-stack seismic dataset, wherein the related seismic facies contain different fluids but similar lithologies than the seismic facies; perform seismic interpretation of the pre-stacked seismic dataset based on the seismic facies and the related seismic facies to generate a digital image of the seismic interpretation; and display the digital image on a user interface.
 17. The computer system of claim 16, wherein the machine learning identifies the seismic facies based on clustering of amplitude variation with offset (AVO) in the pre-stacked seismic dataset, and wherein the seismic facies correspond to AVO clusters.
 18. The computer system of claim 17, wherein the one or more programs, when executed, further cause the computing system to: generate an ensemble of fluid-related AVO responses for the AVO clusters based on a fluid substitution operation; and identify, based on the ensemble of fluid-related AVO responses, related seismic facies across faults indicated in the pre-stack seismic dataset, wherein the related seismic facies contain different fluids but similar lithologies than the seismic facies.
 19. The computer system of claim 16, wherein the pre-stacked seismic dataset includes data relating to a plurality of traces recorded via a plurality of seismic sensors.
 20. The computer system of claim 16, wherein the data relating to the plurality of traces has been pre-processed using at least one of deghosting, multiple removal, or spectral shaping. 